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Abstract: An estimation of the higher twist contribution to the polarized nonsinglet 
structure function gi is given. Its first moment is generally assumed to be small but 
this could be due to a cancellation of its contribution in different Bjorken-a; ranges. We 
estimate the magnitude of (7*"'^ as a function of x in the framework of the renormalon 

method. The calculated higher-twist-correction to gi turn out to be in fact small for all 

(2) 

X values. The correction to the determination of the nonsinglet matrix element d)qg from 
data for transversely nucleon polarization due to the calculated higher twist effects is of 
order 10%. 
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Recently the E143 collaboration measured the transverse polarized structure function 
for proton and neutron |l| and obtained a first approximate result for the nonsinglet 
twist-3 matrix element d)^s- This analysis relies on the assumption that the higher twist 
contributions to gi~^ are negligible. This experimental result is in strong disagreement 
with predictions from lattice calculations 0, and in marginal agreement with predictions 
from QCD sum rule calculations such that possible sources of error have to be looked 
into very carefully. One such uncertainty in the experimental determination of dj^g is the 
size of higher twist contributions to gi~^ . Theoretically they could be large even if they 
give only a small contribution to the first moment of for example if the latter is due 
to a cancellation of twist-4 contributions from different Bjorken-x regions with different 
signs. 

To know the size of the higher twist contributions is also very important for recent 
attempts to obtain AG(x) from the Q^- dependence of gii^x^Q"^) ^. We present a 
renormalon estimate showing that there are in fact cancelations for the first moment of 
Qi"'^. The higher twist contributions to gi{x,Q'^) are nevertheless small, such that the 
mentioned analysis of experimental data is valid. 

Neglecting higher twist contributions the moments of the nonsinglet polarized structure 
function gi{x,Q'^) can, in the framework of operator product expansion, be identified 
with the twist-2 nonsinglet operator matrix elements Ag-^n (we abbreviate gi{x,Q'^) : = 
glM^) - g-M^)): 

Mg,,n{Q^) = f' dxx^-'g,{x,Q^) 

J 

^<;i,n(/i^) (1) 

Here we introduced the factorization scale ^ and the Wilson-coefficient C„. is the strong 
coupling, defined by := Higher twist operator matrix elements are not written. 

To estimate higher twist contributions to eq. (|I]) we use the renormalon-method ||^. 
The forward Compton-scattering-amplitude is calculated in the Borel-plane with an ef- 
fective gluon propagator, taking only one gluon exchange into account. In the Landau 
gauge the effective gluon propagator resums arbitrary many quark loop insertions in the 
gluon line, which is an exact procedure in QED. The gluon propagator reads in the Borel 
representation 0: 

ByaAa.Df^{k)]{u) = 6^" '^^ 1^) , (2) 




C corrects for the renormahzation scheme dependence (C = — | for MS-scheme) and u is 
the Borel transformation parameter. In QCD the restriction to one gluon exchange is an 
exact procedure in the large A'^p-limit [§, |^, [T^, where Np is the number of quark-fiavors. 
The next-to-leading A^i?-order terms are approximated by naive-nonabelianization (NNA) 
|ll|, 0, [T^ , which means to replace the one loop beta-function of QED by the QCD beta- 



function Po = 11 — and corresponds to the replacement Np Np — y. The quality 
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of this approximation has to be checked by comparing the NNA-perturbative-coefficients 
with the known exact ones. 

The IR-renormalon-poles lead to an factorial growth of the perturbative series, which 
has to be interpreted as an asymptotic expansion. In agreement with the general theory 
of asymptotic series it has to be truncated after the minimal term, which determines the 
best accuracy which can be achieved using perturbative expansion. As its dependence 
is power-like, it has been suggested to use it as an estimate of higher-twist contributions. 
Despite the fact that the conceptual basis for this approach is controversial fl^, the pro- 
cedure has given previously reasonable estimates |]l5|, 16, 17 1 and is much simpler than 



genuine higher twist estimates due to lattice or QCD sum rules calculations. 

Finally, let us mention that there are two popular methods of computing the renormalon 
ambiguity in the QCD perturbative series. Apart from the Borel-transformed propagator 
(H) one can use a gluon propagator with a small gluon mass A^. The coefficient of the n-th 
renormalon pole is the same as the coefficient of the A^"'/n(A^)-term in the small gluon- 
mass approach. Despite seemingly different physical meaning, both methods are exactly 
equivalent 0, |TB[. 

The paper is organized as follows: first the Wilson-coefficient for all moments of gi are 
calculated in the NNA approximation using the Borel-transformed gluon-propagator. Then 
the perturbative coefficients are calculated in the NNA approximation and are compared 
to known exact values. Next the IR renormalon ambiguity is calculated and used for an 
estimate of the twist-4 corrections to the polarized structure function. As the calculation is 
done for all moments, the twist-4 corrections are given as a function of Bjorken-x. Finally 
we give a phenomenological analysis of the twist-4 contribution. 

The forward Compton scattering amplitude corresponding to the polarized structure 
function gi reads after an expansion in = 1/x: 
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We have computed the one-gluon-exchange diagrams contributing to deep inelastic lepton- 
nucleon scattering using the effective gluon propagator (Q). 

The result in the Borel plane is expanded in u; and compared to eq. (Q) to read off the 
Borel transformed Wilson coefficient. Note that the Borel transformation of eq. (|^) is given 
by a similar expression with the Borel transformed Wilson coefficient on the right hand 
side. The resulting antisymmetric part is (we use the notation of Muta |l^ and s = Pqu): 
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For the first moment this result coincides with the one obtained by 0. There are IR- 
renormalons appearing at s = 1, 2 and UV-renormalons depending on the considered mo- 
ment of g\. The fact that there are only 2 IR-renormalons appearing is general for the case 
of deep inelastic lepton-nucleon scattering, but not imperative for other processes. The 
1/s-term has still to be renormalized by a corresponding counterterm, which, to lowest 
order, is given by the coefficient in front of the 1/s singularity: 
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where 7^ is the anomalous dimension corresponding to the parton to parton branching 
with helicity transfer. This expression vanishes for n = 1, so that there exists no 1/s 
singularity in the case of the Bjorken sum rule. The factor 2 is due to the fact that we 
defined = agjiAm) instead of as/(27r). The relative sign can be understood as follows: 
The coefficient of 1/s has to be identified with the coefficient of ln(/i^) while the 
anomalous dimension appears in front of ln(Q^//i^) = In(Q^) — ln(/i^). This gives rise to 
the minus sign. 



The Bjorken sum rule has been determined perturbatively up to high orders in [ pl|] 
Writing the result in form of an expansion in Np one gets: 



/ dxgilx^Q"^ 
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where the Riemann zeta functions are given in decimal form and in the last line Np = 3 
was inserted. 

The perturbative corrections to the Bjorken sum rule can be reconstructed from the 
Borel transformed Wilson coefficient Cg^^i in eq. in the NNA approximation by succes- 
sive differentiation with respect to s. Let m be the order of the perturbative correction 
and -B^"*^ the corresponding coefficient. Then 
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This equation remains valid for higher moments if the renormahzed Wilson coefficient is 
inserted on the right hand side. In the MS scheme we get for the NNA approximated 
perturbative coefficients: 



dx gi{x, 
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(8) 



where we inserted Np = 3 in the last line. So the contributions from leading order in Np 
are reproduced exactly, while the nonleading terms are found to be approximated with an 
relative error of about 30% (see also [p!0| ). 

As already mentioned the terms in perturbative series in QCD have the general be- 
haviour to reach a minimum at some order mo of the coupling and to diverge for orders 
m > niQ. So the series has to be truncated at the order mo to get a meaningful result. 
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The ambiguity of this procedure is just of the order of the smallest term and has the 
form of a power suppressed contribution. In the Borel plane this uncertainty is related to 
an IR-renormalon pole, which prevents the unambiguous reconstruction of the resummed 
series. This ambiguity appears as the imaginary part of the Laplace integral: 
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Using the Borel transformed Wilson coefficient (H) 
ambiguity 

ACf { ri(3 + n) 



we get for the leading IR-renormalon 
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The formally undetermined sign is due to the two possibilities to choose the integration 
path in eq. (p!OD avoiding the renormalon singularity. We checked in the case of the Bjorken 
sum rule, that this result is reproduced up to logarithms by the explicit calculation of the 
minimal term in the asymptotic perturbative series. It is tempting to assume that the 
resulting sign of the higher-twist contribution is identical to that of the minimal term in the 
perturbative series generated by the IR-renormalon. In earlier applications this procedure 
lead to results, which agree with experimental data [0 . Accordingly we suppose that the 
sign in eq. (|n|) should be negative, but strictly speaking it remains unknown. 
In an analogous way the nonleading IR-renormalon contribution is calculated: 
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Note that it is not clear at this point whether this contribution is more important than 
the ttg-correction to the leading power correction or not. 

To use this result for an estimate of power corrections to gi as a function of Bjorken-x, 
eq. ([Tl|) and (|12D are transformed back to momentum space 



C'^'\x) = ±^ \4 + 2x + 4x'- + S\x - 1) - 5(5(x - 1) \ (13) 

" Po [ (1 - a;) 

and 
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with 



C«„ = ( rfxx"-iC;«(x) (15) 



and 1/(1 — x)+ is the singularity subtracted contribution to the integral (p!5|). 

While writing up this contribution we learned of an independent calculation |^ for 
this quantity using the equivalent small gluon mass approach. Both calculations agree 
completely. 

The ratio of the moments M*^^ of the twist-4 part of gi and M^^^ of the experimental 
values for gi (containing all twist) reads 



as 1 
Q2' Q4 



+ (16) 

where we used -Bg'^-'n = 1- In this way gl"''^{x,Q') is expressed by a convolution of 
grix^Q'y. 

Q y \yj 
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where C'^^'^ (y) eq. (0) was inserted. 

Using a fit to the available experimental data for gi |2^, this expression indeed leads 
to very small twist-4 corrections to gi{x) (see fig. 0, |^). This result is not significant for 
the small x region, where the large- A^^p limit is no longer justified. 
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Figure 1: The experimental fit for gi(x, AGeV'^) (full line) [^. The estimate for the twist-4 
correction is given explicitly multiplied by a factor of 10 (dotted line), showing a change in sign 
at X ^ 0.43. The correction was added and subtracted from the experimental values (dashed 
lines), so that the theoretical prediction of the twist-2 part of gi has to lie within this region 

(A^ = 200 MeV, = 4 GeV and Nf = 4). 



From our general result for 5fJ"''^(x, Q^) eq. {\L7\} we can also reproduce the renormalon 



prediction for the nonsinglet twist-4 matrix element |TB 
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in good agreement with our QCD sum rule result of — 0-0094 GeV Another renormalon 

estimate based on Fade approximants [^| gave a value of GeV ^ ^^^^ same way 

as we did in [|1^] we checked the cancellation of the above IR-renormalon contribution by 
the UV-renormalon contribution to the genuine twist-4 correction. Note that in addition 
to the twist-4 power correction to the Bjorken sum rule there exist twist-2 and twist-3 
matrix elements, which are also power suppressed by 1/Q^ p6| . 

Finally we give an estimate for the error made in the determination of the nonsinglet 
twist-3 matrix element c?'-^-' due to higher-twist corrections. The third moment of the 
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Figure 2: The relative magnitude of the predicted twist-4 contribution to gi with respect to 
the experimental fit. Aj^ = 200 MeV, = 4 GeV and Nj = 4. 



nonsinglet structure function g2{x, Q^) := g2{x, Q^) — g2{x, Q"^) is given by 

jjxx^g^ix^Q') = + + (19) 

where a*^^^ is a pure twist-2 nonsinglet operator matrix element. So the twist-3 nonsinglet 
matrix element d*-^-' can be obtained by measuring g2{x,Q'^) and a*-^-', which in turn is 
determined by the third moment of gi{x, Q^): 

£dxx'g,{x,Q') = + (20) 

If one uses a*^^^ from the total measured gi{x, Q'^), higher twist corrections to gi{x, Q"^) are 
not distinguished from the twist-2 part. This leads to an error in the determination of 
(i*-^^ of the order of the higher-twist correction to gi{x^ Q^) Using our estimate for the 
nonsinglet twist-4 contribution to gi{x,Q'^), we obtain the following uncertainty for a*^^^ 

2 dx x^g'^Ux, ^ 0.0003 - Afa^^)) (21) 
Jo 
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The main contribution to the third moment of g2{x, Q"^) comes from a'-^-*, so that the above 
uncertainty leads to an error for rf*^^^ of about 10%, which is substantially smaller than the 
present experimental error bars. We conclude that the mentioned procedures to determine 
S-"^^ and AG(a;) from experimental data are not invalidated by higher-twist effect. 
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